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We carry out a systematic study of a novel type of chaos at onset ("soft-mode turbulence") based on nu- 
merical integration of the simplest one dimensional model. The chaos is characterized by a smooth interplay of 
different spatial scales, with defect generation being unimportant. The Lyapunov exponents are calculated for 
several system sizes for fixed values of the control parameter e. The Lyapunov dimension and the Kolmogorov- 
Sinai entropy are calculated and both shown to exhibit extensive and microextensive scaling. The distribution 
functional is shown to satisfy Gaussian statistics at small wavenumbers and small frequency. 

Spatiotemporal chaos (STC) is a subject of considerable ex- 
perimental and theoretical importance and occurs in a wide 
variety of driven, dissipative systems [|l]-^||. Such chaotic be- 
havior in spatially extended systems is extremely difficult to 
characterize quantitatively, as the dynamics involves a large 
number of degrees of freedom. The most common and use- 
ful tool for the characterization of chaos is given by the Lya- 
punov exponents {A;}. Knowledge of this Lyapunov spec- 
trum permits one to estimate the number of effective degrees 
of freedom of the system (i.e., the dimension of the attractor), 
using for example the Kaplan- Yorke |Q| formula for the Lya- 
punov dimension D(L), where L is the linear system size. It 
also permits one to test the important concept of extensivity 
of chaos, defined as the case in which lirriL-jao D(L) ~ L d , 
where d is the spatial dimension of the system An inter- 
pretation of extensive chaos is that the whole system can then 
be thought of in some sense as the union of almost indepen- 
dent subsystems. This was originally proposed by Ruelle [Q], 
who argued that widely separated subsystems of a turbulent 
system should be weakly correlated, so that the spectrum of 
Lyapunov exponents would be the union of exponents asso- 
ciated with each of the subsystems. The question is closely 
related to the fundamental problem of ergodicity of nonequi- 
librium systems. If the chaos is extensive and each subsystem 
evolves in time practically independently of the others, then 
in a steady (non-transient) chaotic state the time average is 
equivalent to the ensemble average and the system should be 
ergodic. Much work has focused on attempting to character- 
ize spatiotemporal dynamics in these terms (see, e.g., [Q-|ll]]). 
However, in spite of the fundamental importance of the ques- 
tion practically all the results are related just to a few discrete 
coupled map lattices fll2|-|l4|] and two continuous systems, 
namely the complex Ginzburg-Landau (CGL) and Kuramoto- 
Sivashinsky (KS) equations (see, e.g., refs. Jl3| , p3p^ ]). This 
is partly due to the computational complexity of the problem 
but primarily to the lack of simple models exhibiting the req- 
uisite chaotic behavior. It would therefore be of considerable 
interest to characterize quantitatively other types of STC. 

Recently attention was drawn to the existence of a new wide 



class of systems displaying such a behavior [|17|-|2C|]. Their 
properties are qualitatively different from those of the CGL 
and KS models. In contrast to both these models the chaos 
is associated with smooth, random long-wavelength modula- 
tions of a short-wavelength pattern, with defect generation be- 
ing unimportant. The short-wavelength pattern arises due to 
a single supercritical bifurcation of the Turing type such as 
occurs in Rayleigh-Benard convection. The long-wavelength 
modes belong to a Goldstone branch of the spectrum origi- 
nated in a broken continuous symmetry. The symmetry makes 
the system degenerate to the extent that instead of a single, 
unique spatially uniform state, it has a continuous family of 
equivalent spatially uniform states, which may be obtained 
from each other by the symmetry transformation. This sym- 
metry, which is additional to the trivial groups of translations 
and rotations, can be one of many different types. For this 
reason the STC in question is quite a common phenomenon 
and occurs, for example, in electroconvection in liquid crys- 
tals [ [l7||l8| ], in convection in a fluid with stress-free bound- 
ary conditions [{2l|-|24}|, etc., see ref. $ldj\ for further discus- 
sion. The chaos observed in such cases may be interpreted as 
a macroscopic dynamical analog of second order phase tran- 
sitions, where the order parameter is related to the amplitudes 
of turbulent modes. Due to this analogy it has been called soft- 
mode turbulence (SMT) jl8[]. The simplest model exhibiting 
SMT was introduced by Nikolaevskii [ p5| , p6| ] to describe lon- 
gitudinal seismic waves in viscoelastic media. In what follows 
we exploit the simplicity of this model to shed light on general 
features of this new type of STC. 

We present a detailed systematic study of the Nikolaevskii 
model, including calculation of the Lyapunov exponents, frac- 
tal dimension, etc. We show that the Lyapunov dimension and 
Kolmogorov-Sinai entropy are extensive quantities, which 
supports the validity of the Ruelle's ideas for SMT. We also 
show that the power spectrum can be described by Gaussian 
statistics for small wavenumbers and frequencies, in agree- 
ment with a general argument of Hohenberg and Shraiman 
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The model is defined by the following partial differential 
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equation for the real scalar field v(x, t) (longitudinal mode of 
the displacement velocity in the original formulation) [E5pq]: 



dv 
dt 



dx 2 
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with < x < L and periodic boundary conditions. This 
model has two control parameters, e (the distance from on- 
set) and L, in contrast to, e.g., the KS model, where the only 
non-trivial control parameter is the system size L. An essen- 
tial feature of the model is that even at small e it cannot be 
reduced to any e-free form [20,27j], which makes the hierar- 
chy of characteristic scales much more complex than those 
in the KS and CGL models Eq. (1) may be re- 

garded as a generalized Burgers equation and shares with it 
the same group of symmetry, namely trivial symmetries un- 
der shifts of the spatiotemporal coordinate system, and non- 
trivial invariance with respect to the Galilean transformation 
v(x,t) — > v(x—v t,t)+v , where v is an arbitrary constant. 
The Galilean invariance plays the role of the above specified 
additional symmetry, generating for Eq. (1) the continuous 
family of solutions v = v . The Nikolaevskii equation admits 
a continuous set of spatially periodic, stationary solutions, 
whose instability has been proved analytically [^8j. Computer 
simulations [ |l9| ] of an equivalent version of this model for an 
order parameter u(x,t) = J v(x,t)dx showed that even at 
extremely small e = 10~ 4 , the system exhibits STC. How- 
ever, this simulation did not provide any quantitative results 
about the STC. The only result of this kind is in ref. [|3C|], in 
which just a single quantitative characteristic, namely the de- 
pendence of the mean amplitude of chaotic patterns on the 
control parameter, was studied. In this Letter we span the gap 
in our knowledge of this new type of STC, providing a de- 
tailed quantitative description of its most important properties 
based on numerical integration of Eq. (1). The simulations 
were carried out using the pseudo-spectral method combined 
with a fourth-order predictor-corrector integrator, for several 
different values of L and two values of e (0.2 and 0.5, re- 
spectively) [pTj]. The Lyapunov exponents were calculated 
by linearizing the equation along the trajectory, performing a 
re-orthonormalization after a few integration steps to prevent 
the largest Lyapunov exponent from swamping all the others 
J32I]. A typical pattern v(x, t) as a function of space and time 
in the steady chaotic regime is shown in Fig. 1. The time 
averaged power spectrum (\v(k)\ 2 ) obtained over a time pe- 
riod of T = 10 4 for system size L = 78 with e = 0.2 is 
shown in Fig. 2. As can be seen, the dominant modes occur 
in the vicinity of k = ±1. Note, however, the smaller peak 
near zero wavenumber, which arises from coupling between 
unstable short-wavelength modes centered about k = ±1 and 
the slowly decaying modes from the Goldstone branch of the 
spectrum centered at k = (using terminology based upon 
the linear stability analysis of the spatially uniform solution). 
The power spectrum for u(k) = v(k)/ ik (shown in the insert 
in Fig. 2) has the dominant peak near k — 0, as was origi- 
nally demonstrated in [|l9||. One would also expect on gen- 
eral grounds [pi] each Fourier transform variable v(k,uj) to 



be governed by a Gaussian probability distribution functional, 
exp(— D(\v(k,u})\ 2 ), for small k andw. We have verified this 
for several values of k and u>. In particular, for uo — 0.6 the 
Gaussian distribution holds for < k < 2 (see Fig. 3). 
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FIG. 1. Typical space-time configuration for v(x, t) in the steady 
chaotic regime is shown here. The configuration shown has evolved 
from random initial condition within a system size L — 78 and 
e = 0.5. 




FIG. 2. Time average power spectra {\i>(k)\ 2 } and {|ii(fc)| 2 ) in 
fc-space are plotted. The system size is L — 78 and control parameter 
e = 0.2. Note that the power spectrum has symmetry with respect to 

k-t-k. 



As has already been emphasized, the key question for any 
system exhibiting STC is whether it can be represented as a 
union of weakly correlated subsystems. If this is the case, 
then the spectrum of Lyapunov exponents for the entire sys- 
tem should be intensive in the sense that A, is a function only 
of the intensive index i/V, i.e., A; = f(i/V), where V stands 
for the volume of the system [H]. The question is not trivial 
for the type of STC considered here, because of the impor- 
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tance of the long-wavelength modes and the divergence of the 
two point correlation length £2 p3| ] as e — > [p0|], To an- 
swer this question, a detailed study of the Lyapunov spectrum 
for the Nikolaevskii model was conducted. The results are 
shown in Fig. 4, where the number of Lyapunov exponents 
greater than a particular value A;, scaled by the system size L, 
is plotted versus A; for e = 0.5. A similar curve is found at 
e = 0.2, but the maximum positive eigenvalue is now smaller 
(one expects this eigenvalue to vanish as e — > 0.) The inten- 
sive nature of the Lyapunov density is evident. In this case, 
one also expects the fractal dimension D(L) of an attractor 
to be extensive for large enough L (extensive chaos), as was 
first shown by Manneville [15| for chaotic solutions of the 
Kuramoto-Sivashinsky equation. We have checked this for 
the Nikolaevskii model using the Kaplan- Yorke formula [|j] 
for the Lyapunov dimension, 



K 



D(L) = K + J2 x i/\ X 



fc+H 



(2) 



where the integer K is the largest integer such that the sum 
of the first K Lyapunov exponents is nonnegative. We have 
also calculated the Kolmogorov-Sinai entropy H (L) from the 
definition 



ff(L) = 5> 



(3) 



where the sum is over the positive Lyapunov exponents. The 
Kolmogorov-Sinai entropy tf32[ ] is a measure of the mean rate 
of information production in a system, or the mean rate of 
growth of uncertainty in a system subjected to small perturba- 
tions. We find that for large enough L both D(L) and H(L) 
are extensive. Our results for the Lyapunov dimension D(L) 
are shown in Figure 5 for e = 0.5. The same behavior is found 
at e = 0.2, but with a different slope (naturally the slope of the 
straight line D(L) at e = 0.2 is smaller than that at e = 0.5). 
In addition, we find that the upper index i + in (3), which cor- 
responds to the smallest positive Lyapunov exponent, is also 
proportional to L. 
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FIG. 3. Probability density distribution of power spectrum v(k) 
for k — 1.5 with system size L = 78 and e = 0.5. Notice that 
the validity of the Gaussian distribution for v(k) (see the main text) 
implies an exponential distribution f(z) — De~ Dz for the variable 
z = \v(k)\ 2 . In this case we find the value D ~ 47. 

The important characteristic of STC is the dimension cor- 
relation length £5 This length is defined as £g = <5 _1 / d , 
where S = Urn l _> ^ D (L) / L d . It can be thought of as the 
"radius" of a volume that contains one degree of freedom, or, 
as the linear size of the subsystem described above. The value 
of this dimension correlation length for Eq. (1) is £5 = 3.0 for 
e = 0.5 and £5 = 3.3 for e = 0.2. In contrast to £5 the above- 
mentioned two point correlation length £2 governs the spatial 
decay of the correlation function. In general these two lengths 
are different [[3j]. We found that £ 2 — 4.9 for the Nikolaevskii 
model at e = 0.5 and £2 — 5.6 at e = 0.2. 

Finally, Tajima and Greenside [[34]] have recently found for 
the one dimensional Kuramoto-Sivashinsky model that D(L) 
is also "microextensive." Namely, they found that if one in- 
creases L by a small amount SL, with SL <C £5 (SL = 0.8 
in our simulation), one finds that D(L) satisfies the same lin- 
ear relationship as that characterizing extensive chaos. We 
have examined this for two different domains of L for our 
model and found that microextensivity holds for both D(L) 
and H(L). Our results for microextensivity for D{L) are 
shown in the insert in Fig. 5. 




FIG. 4. Here iV>i is the number of Lyapunov exponents greater 
than a particular value \i . We plot N>i / L (scaled by the system size 
L) vs Ai in the case e = 0.5. 

In conclusion, our detailed study of this new type of STC 
based on the numerical integration of the Nikolaevskii model 
shows that for sufficiently large system size the chaos is both 
extensive and microextensive. We also found that the system 
satisfies Gaussian statistics at sufficiently small wavenumbers 
and frequencies. We believe these results are quite general 
and reflect intrinsic features of this type of STC, rather than 
specific peculiarities of the model. 

There are several interesting questions to investigate in the 
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limit e — > 0, including the dependence of quantities such as 
the correlation lengths and Lyapunov exponents on e as well 
as the possible scaling of the power spectrum. This study is in 
progress and will be reported elsewhere. 
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FIG. 5. Lyapunov dimension D(L) vs system size L for e = 0.5. 
SL = 15.5 for extensive, and 8L — 0.8 for microextensive. 
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